Bose-Einstein Condensates in Strongly Disordered Traps 
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A Bose-Einstein condensate in an external potential consisting of a superposition of a harmonic 
and a random potential is considered theoretically. From a semi-quantitative analysis we find the 
size, shape and excitation energy as a function of the disorder strength. For positive scattering 
length and sufficiently strong disorder the condensate decays into fragments each of the size of the 
Larkin length C. This state is stable over a large range of particle numbers. The frequency of the 
breathing mode scales as Ij . For negative scattering length a condensate of size £ may exist as 
a metastable state. These findings are generalized to anisotropic traps. 

PACS numbers: 03.75Hh, 03.75Kk 



Bose-Einstein condensation (BEC) is of great inter- 
est in a wide variety of systems including superfluidity 
and superconductivity. Its perhaps cleanest realization 
is found in dilute ultracold gases of alkali atoms conflned 
in magnetic or optical traps which have been studied in 
great detail both experimentally [l|, Q and theoretically 
H, |3| (and references therein). 

More recently these investigations have been extended 
to traps formed by a superposition of a harmonic and a 
random potential [1, @, 01 • Whereas for harmonic traps 
the behavior of the condensate, e.g. its size, shape, ele- 
mentary excitations etc. is quite well understood, this 
is not the case for traps including a finite amount of 
disorder. Most of the theoretical investigations of BEC 
in random potentials use as a starting point a descrip- 
tion in momentum space (e.g. Bogoliubov transforma- 
tion or Beliaev-Popov perturbation theory) where the 
study of the typ ical disorder effects is notoriously diffi- 
cult 

EBEEf- Disorder is there treated perturbatively 
with as a small parameter, where ^ and L denote the 
healing and the Larkin length, respectively. As we will 
show below in the case of sufficiently strong disorder this 
parameter is of order one and perturbation theory breaks 
down (but see [l^). 

In the present article we present a semi-quantitative 
analysis of BEC on the level of Larkin-Imry-Ma argu- 
ments 131 llJ] to determine the size, shape and the fre- 



quencies of the breathing mode in weakly interacting 
Bose-Einstein condensates in disordered traps. In addi- 
tion to the oscillator and the scattering length the Larkin 
length appears as a new relevant length scale which dom- 
inates the properties of the condensate for weak interac- 
tion and strong disorder. Strong disorder reduces the 
anisotropy of the condensate. For attractive interaction 
the condensate may exist in a metastable state. 

The model. Starting point is the Hamiltonian of inter- 



acting bosons 



2m 



C/(x) 



^-t*)*. (1) 



Here ^'^^(x) and ^'(x) are the creation and annihilation 
operators of the Bose field with J d'^xj^'p — N for the 
particle number, m denotes the mass of the bosons and a 
the scattering length. The external potential represents 
a superposition of the harmonic trap and a Gaussian ran- 
dom potential 



C/(x) = |^c^,2x2 + [/di3(x). 



(2) 



i=l 



The random potential is characterized by 

(t/dis) = 0, ([/di.(x)C/dis(x')) - k'<5(x - x'). (3) 

With this ansatz we assume that the correlation length 
of the disorder is smaller than all other length scales. We 
will briefly discuss other cases at the end of this article. 

In what follows we consider the state of the system at 
zero temperature. We will first examine the case of an 
isotropic 3-dimensional trap with uji ^ lu2 ~ lo^ — luq. 
Assuming a compact condensate cloud of the radius R 
the condensate energy per particle e can be written as 



£(i?) 
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(4) 



Here we introduced the oscillator length (the size of the 
harmonic oscillator ground state) i = {h/ imujo))^/'^ . C 



13 



denotes the Larkin length C ~ IGttH'^ / (27m'^ k^) 
14l | in d = 3 dimensions. N is the total particle number. 
In the following we will set h = m = 1. In this notation 
the healing length in the absence of disorder reads ^ = 

The first and the second term in ([4]) describe the kinetic 
energy and the interaction of the particles, respectively. 
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Both terms favor the spreading of the condensate pro- 
vided the scattering length is positive. If the radius of 
the cloud is sufficiently large, namely 



i? > i?a = 3a7V, 



(5) 



the interaction is negligible in comparison with the ki- 
netic energy. Neglecting the inhomogencity of the con- 
densate, condition ([5]) can be rewritten as ^ ^ i?. The 
third and the fourth term describe the oscillator potential 
and the typical potential well from the disorder, respec- 
tively, both cause its confinement. In the last term we 
have taken into account that the typical value of the fluc- 
tuation of the potential in a region of linear size R scales 
as kR^^/^. This term is larger than the kinetic energy 
only on scales R > C. Similarly to the interaction, at 
sufficiently large radius of the cloud. 



R:$> Rc^ei£/c)^/\ 



(6) 



the disorder can be neglected in comparison to the har- 
monic potential. 

In general the center of an attractive domain formed by 
static fluctuations of the random potential may not coin- 
cide with the center of the harmonic trap. Nevertheless, 
our estimates are correct since it is either the harmonic 
trap or the disorder which leads to the localization of 
the condensate. Experimentally localization by disorder 
is often observed by a sudden decrease of the oscillator 
frequency. In this case the condensate will be localized 
in the coarse grained potential well closest to the origin. 

Condensate size. Below we determine the equilibrium 
size i?o of the condensate as a function of the particle 
number and the disorder strength from the equilibrium 
condition de/dR = 0. In the case of weak disorder and 
relatively small particle number, such that £ <g; £ and 
N <C £/(3a), we can ignore both the disorder and the 
interaction and hence, from ^ Rq k, £. At large particle 
number TV > ^/(3a), the effect of the interaction sets in 



and i?o approaches the Thomas- Fermi radius [15 1 



(7) 



More interesting is the case of strong disorder, such 
that the oscillator length is much larger than the Larkin 
length, £ ^ C. Then, for small particle number, the com- 
petition between the kinetic energy and the effective po- 
tential well resulting from the random potential on scale 
R gives 



Ro^ C, Na-^C. 



(8) 



Thus the size of the condensate is equal to the Larkin 
length. We estimate the healing length as ^ '--^ 
■\/ C^/Na ^ C from which it is plausible to conclude 
that the condensate is strongly depleted and the super- 
ffuid stiffness vanishes [sl. 



For larger particle numbers the interaction becomes 
important and the minimum of (|4]) is given by 



/ 9 \ 2/3 , / £\ 

Ro{N)^i^-Naj £1/3, /:«aiV«f[-j 
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(9) 



The (average) healing length is estimated from ^ ~ 
^Rl/Na - {C/Naf/^Ro < i?o from which we con- 
cludes that the condensate is weakly depleted only and 
the superfluid stiffness remains flnite [S!]. At the up- 
per boundary for N the radius reaches Rc- Finally, for 
even larger particle number the radius crosses over to 
the Thomas-Fermi radius. The crudeness of our energy 
expression ([4]) does not allow accurate determination of 
the numerical pre-factors for all quantities derived from 
(|4]). If the disorder decreases, the Larkin- length C in- 
creases and the interval where ^ applies diminishes and 
eventually vanishes when C reaches the oscillator length 
£. In this case only the cross-over at Na ~ £ discussed 
previously survives. 

Fragmentation. Next we consider the instability of a 
compact condensate against fragmentation. For N <C 
C/a the ground state energy for increasing N decreases 
as E K, —N/iQC"^) until it reaches a plateau E k, Eq at 
N>No~ C/a where Eq « -2/{27Ca) (compare ®). 
Thus, at > A^o, the division of the condensate into 
N/Nq fragments is energetically favorable. The energy 
of the fragmented condensate is « {N/No)Eo{No) w 
—NjC? < Eo{N). Thus the fragmented condensate cor- 
responds to the true ground state. The total volume 
of all fragments scales like ~ NaC^ and is hence by a 
factor {C/{Na)) smaller than the volume of the com- 
pact state. The healing length ^ in the fragments is 
i = {C'^/{ma)f/'^ where N < Nq, and hence ^ > C. Ac- 
cording to 01 , in the region where ^ > £ a finite fraction 
of the bosons are not in the condensate. The superfluid 
stiffness of the fragmented state most likely vanishes. To 




FIG. 1: State of dense fragments. Each fragment includes 
£/a particles. The minimal condensate radius is Rf ~ 

(iVa)i/^£2/^ 

describe the cross-over to the Thomas- Fermi behavior we 
consider a state of dense fragments with the total radius 
Rf ^ {NaY^^ L^^^ . Rp denotes the minimal size of of 
the fragmented state. If there is no confining parabolic 
potential the fragments may be distributed over a much 
larger area to take advantage of the tails of the Gaus- 
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FIG. 2: Condensate radius as a function of the particle num- 
ber. The bold hne corresponds to the equilibrium fragmented 
state, the dashed line corresponds to the metastable compact 
state. 



sian distribution of disorder. In a fragmented state with 
radius Rp the oscillator energy per particle is of the or- 
der Rp£~'^ which has to be compared with the energy 
per particle —C~^ from the disorder. This gives for the 
cross-over particle number 



Nr.. 



£{£/cy 



(10) 



This number is much larger then the particle cross-over 
number between ([9]) and ([7]) provided £ > C. 

The compact state described by ^ may occur as a 
metastable state which lives for a short period of time 
after a confining harmonic trap is switched off (or the os- 
cillator frequency lOq is strongly reduced) . On larger time 
scales the condensate will lower its energy by decaying 
into ^ Na/C fragments each of the size C. This decay 
into fragments may explain the striped phases found in 
cigar-like traps Q . 

Breathing mode. Next we consider the breathing mode 
frequency of the radial oscillations of the condensate 
around the equilibrium configuration from 



mul = e"{Ro). 



(11) 



Without disorder we find from ((4]) and (fTT|) ujb — 2a;o 
in the non-interacting and uj = y/buJo in the interacting 
case, respectively, in agreement with previous results [l6| . 
In the disordered case we obtain for Na <C £ <C ^ 



UJB 



1 
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UJo > t^o 



(12) 



and for C <C aN <C £{£/C)^^'^ under the assumption of a 
compact non-fragmented condensate 



UJB 



UJo » Wo. 



(13) 



Contrary to the pure case the interaction changes the 
radial oscillations of the condensate strongly. 

In the true equilibrium fragmented state the breathing 
mode will be controlled by the weak interaction between 
the fragments. In this case the frequency gap will prob- 
ably decrease as ~ Rp^- 



Negative scattering length. So far we assumed that the 
scattering length a is positive. In the case of negative 



scattering length as e.g. in the case of ^Li atoms 
there is only a metastable state of finite radius Rc follow- 
ing from de/dR\R=R^ = 0. This state becomes unstable 
at a critical particle number Nc. The condition for the 
disappearance of the metastable minimum follows from 



de{R) _ dhjR) 
OR ~ ~dW 







R=Ra 



(14) 



In the pure case £ ^ £, condition gives for the 
critical particle number N^ « 0.12^/|a| in agreement with 



The strong disorder decreases the radius of the 
metastable droplet as it follows from the equilibrium con- 
dition 



1/2 



4 N 



3/2 



(15) 



Here Nc denotes the critical particle number beyond 
which the metastable state disappears 



8 £ C, 
Nc = —n ~ 0.033--. 

243 a a 



(16) 



The radius of the metastable droplet in the region N < 
Nc decreases monotonically from Rc ^ C for N <^ Nc to 
Rc ~ §£ for N w Nc. 

Anisotropic traps. So far we assumed that the trap is 
isotropic whereas experimentally often anisotropic traps 
are considered. Our results can be easily extended to 
these cases. For brevity we will consider here only the 
strongly anisotropic situations (i) uji ^ uj2 = uJq uj^ — 
uj± [d = 2) and (ii) uji = uJq uj2 — uj^ = cjj_ {d — ^) 
corresponding to a pancake and a cigar-shape like trap, 
respectively. We assume that the the oscillator length 
obeys £±_ — (/i/mcjj^)^/^ ^ ^, £ so that the extension of 
the condensate in the transverse direction is given by £±_ . 
The extension i?o in the longitudinal direction follows 
then from the minimum of the energy per particle 



1 

2R^ 



ZaN 



2R'^£^, 



R^ 

-d + 2^4 



4-d> 



1/2 ' 



d= 1,2. 
(17) 

Here we have introduced the d— dimensonal Larkin length 
Cd = (££^~'*)^/(^^'*) and assumed that the correlation 
length b of the random potential is small compared with 
the transverse width of the trap, b <^ £±. In the opposite 
case b ^ £±, £± has to be replaced by £']_/b in the defi- 
nition of Cd- From pT]) we can read off the generalized 
cross-over radii 



l/(2-d) 



(18) 



For R^ Ra the interaction and for R ^ Rc the disorder 
can be neglected, respectively. For a d- dimensional trap 
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we find in analogy to ([8]) for small particle numbers 

Ro^Cd, Na< e± ie±/Cdf-^ . (19) 

In this region the healing length £, > Cd such that a 
finite fraction of particles is not in the condensate and 
superfluidity is suppressed. For larger TV Q has to be 
replaced by 

Ro « £^ {aN/i^f" {Cd/t^)^^-"'^''' . (20) 
This formula is valid as long as 

„ , (4-d)(d+2) 

aN < ii_ {l^/lf-'' {l/Cd) (21) 

For even larger N one find the generalized Thomas- Fermi 
behavior 

Rtf{N) « li_ (aiV^V^i) ^ (22) 

However, we have to consider here again the instability 
of the compact domain state described by (|^D|) against 
fragmentation. The energy of the compact state scales 
as — (l/a^j^)(£j^/£rf)''~'' whereas that of the fragmented 
state scales as —N/C\. Thus the fragmented state has 
the lower energy until the radius of the set of dense frag- 
ments 

RF{N)^i^{NaCytiYl'' (23) 

reaches a value such that ~ R'p/i^ which happens 
at 

Nco ~ (^i-^^'-'^a) i£/Cf+'' (24) 

which generalizes pUj) . but is more realistic, especially at 
d=l. 

In conclusion, we demonstrated that the disorder 
strongly influences the size, shape and structure of the 
condensate cloud. Sufficiently strong disorder leads to 
the localization of the condensate fragments of linear size 
of the Larkin length. This result agrees with the conclu- 
sions about the suppression of the transport in elongated 
Bose-condensates derived in The increasing number 
of particles enhances interaction and may delocalize the 
condensate. Such a delocalization transition is especially 
realistic in a strongly prolonged cigar-shape condensate. 
If the atoms in condensate attract each other, the dis- 
order decreases the critical number of particles at which 
collapse is developed. 
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